The Effects of Symmetries on Quantum Fidelity Decay 
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We explore the effect of a system's symmetries on fidelity decay behavior. Chaos-like exponential 
fidelity decay behavior occurs in non-chaotic systems when the system possesses symmetries and 
the applied perturbation is not tied to a classical parameter. Similar systems without symmetries 
exhibit faster-than-exponential decay under the same type of perturbation. This counter-intuitive 
result, that extra symmetries cause the system to behave in a chaotic fashion, may have important 
ramifications for quantum error correction. 
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FideHty decay as a possible signature of quantum chaos 
was introduced by Peres while exploring irreversibil- 
ity in quantum mechanical systems. Though the overlap 
between two initial states undergoing equivalent evolu- 
tion remains constant with time, it decreases with time 
if the Hamiltonian effecting one of the systems is sightly 
perturbed. The behavior of the decrease in overlap, or fi- 
delity, depends on whether the evolution is the quantum 
analog of a chaotic or non-chaotic classical system. 

Further explorations have distinguished realms of fi- 
delity decay behavior based on the chaoticity of the cor- 
responding classical system, perturbation strength and 
type, and the initial system state. Weak perturbations of 
the Hamiltonian, such that perturbation theory is valid, 
exhibit a Gaussian fidelity decay Stronger pertur- 

bations, in the Fermi golden rule (FGR) regime, exhibit 
exponential fidelity decay with a rate determined by the 
strength and type of perturbation The rate of 

this exponential generally increases as the square of the 
perturbation strength j3] and may saturate at the under- 
lying classical systems' Lyapunov exponent 0, @1 or the 
bandwidth of the Hamiltonian 3] . 

The fidelity decay of classical-like coherent states for 
quantum analogs of non-chaotic systems 0, 0, 0| can be 
Gaussian ..llj , faster than the chaotic system exponen- 
tial, or power-law ^2 depending on the initial state Q 
and the effect of the perturbation on the classical orbits 
on which the state is centered jlSl . The fidelity decay 
of systems states and perturbation strengths 
[r^ at regime edges have also been explored. 

The quantum fidelity decay of an initial state is 



(1) 



F{t) = U,\U-%\^.)\^ 

where U is the unperturbed evolution. Up — C/e~"^^ is 
the perturbed evolution, e is the perturbation strength. 



and V is the perturbation Hamiltonian. For chaotic sys- 
tems, random matrix theory (RMT) gives the FGR ex- 
ponential fidelity decay rate F as a function of e and the 
perturbation Hamiltonian eigenvalues Aj 



F = e^i^. 



(2) 
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where = N-^ Af and N is the Hilbert space di- 
mension. 

Fidelity decay is related by a Fourier transform to 
the local density of states (LDOS) H [3 77(A0) ^ 
\{vm\v'n)\^ , where = (j)m — (f)'ni is the difference be- 
tween unperturbed and perturbed eigenangles given by 
the eigenvalue equations: U\vm) = exp(— i(^m)|'ym) and 
UpWn) = 6xp(— The LDOS provides a measure 
of how local is the perturbation. For complex systems 
the LDOS is typically Lorentzian |l9j |. the perturbation 
transfers probability to far reaches of the system basis. 
The Lorentzian width F gives the exponential chaotic fi- 
delity decay rate. 

Recently Rossini, Benenti, and Casati (RBC) [l^l re- 
ported numerical evidence of exponential fidelity decay 
in non-chaotic systems due to 'quantum' perturbations 
(though the rate may not be that expected for chaotic 
systems) . Quantum perturbations are not tied to a clas- 
sical system parameter and are applied multiple times 
during each map iteration, i.e. after each basic gate in 
the simulation of the dynamics on a quantum computer. 
The cause of this behavior, RBC explain, is the non- 
locality of such errors which allow for direct transfer of 
probability over a large distance in phase space. Mean- 
ing, the LDOS has significant amplitude even for large 
A</). This result is of particular importance for suggested 
exploitations of the slower exponential decay to stabi- 
lize quantum computation [2lLl2^. and for fidelity decay 
studies on a quantum computer .20. . .23.. i24 | . 

In this paper we numerically demonstrate that pertur- 
bations not tied to a classical parameter of non-chaotic 
systems can exhibit exponential decay even when applied 
after every map iteration, as generally done in fidelity 
decay studies. We choose this perturbation application 
scheme because 1) there is no unique gate sequence for 
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implementing an operator 2) some of the operators we 
discuss cannot be implemented efficiently on a quantum 
computer 3) perhaps most importantly, when the pertur- 
bation is applied at numerous points during the map it- 
eration the composite effect on the map is likely random. 
Thus, the exponential decay is due to the basis of the 
perturbation as found in 0. By applying the perturba- 
tion every map iteration we show that the decay behav- 
ior is due to a different phenomenon: symmetries in the 
system Hamiltonian. When no symmetries are present 
the perturbation is local, the LDOS is Gaussian, and the 
fidelity decay is faster than exponential even when the 
perturbation is not tied to a classical system parame- 
ter. We emphasize the counter-intuitive nature of this 
claim, namely, by adding symmetry, the system's fidelity 
behavior behaves in a more chaotic fashion. 

Our numerical study centers around the quantum 
kicked top (QKT), a system used in many studies of 
quantum chaos in general "25] and fidelity decay in par- 
ticular 0,0,110. The QKT 26] Floquet operator is 

Uqkt = e^"^''/2g-^fcJ,?/2j^ 

where J is the angular momentum of the top and J are 
irreducible angular momentum operators. The Hilbert 
space dimension of the top is iV = 2 J -f 1 and the rep- 
resentation is such that is diagonal. The chaoticity 
of the QKT depends on the kick strength, k. The QKT 
is non-chaotic for k ^ 2.7, has chaotic and non-chaotic 
regions for 2.7 ^ fc < 4.2, and is fully chaotic for k > 4.2 
y. The QKT has different symmetry sectors based on 
its angular momentum. For all J the QKT has conserved 
parity with respect to 180° rotations about y. For even J 
the subspace even with respect to y has conserved parity 
with respect to 180° rotations about x. 

In his original fidelity decay studies, Peres was careful 
to account for the symmetries of the QKT 0] by block 
diagonalizing and performing simulations using only one 
block. While some later authors have done the same 0, 
others 0, Q] have used the complete QKT without ap- 
parent discrepancies. Here we show that the presence 
of symmetries in a quantum system can dramatically 
change the fidelity decay behavior. 

We choose a 'quantum' perturbation relevant to quan- 
tum control studies 

f/p = n;iiexp(-ieai/2), (4) 

where Uq is the number of qubits and ct^ is the Pauli spin 
matrix. This corresponds to a collective qubit rotation 
about the z-axis by angle e and is a model of coherent 
far-field errors 27]. We apply the perturbation only after 
a complete map iteration. 

Figure shows the fidelity decay of the QKT for dif- 
ferent values of k with the above perturbation, e = .2, 
and random initial states. Values of k are shown for reg- 
ular, mixed, and chaotic QKT but the fidelity decay is 
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FIG. 1: (Color online) Fidelity decay for QKT with J = 
255.5, and A; = 1 (o), 3 (□), 5 (o), and 7 (A) for the collec- 
tive bit z rotation with e = .2. Even the non-chaotic QKT 
decays are exponential though they do not follow the random 
matrix theory rate given in Eq. |5| (dashed line) . The inset 
shows the fidelity decay for k — 1 (x) and 10 (o) for pertur- 
bation strengths e — .2, .3, and .6 (top to bottom). As the 
perturbation strengthens, the invariant subspaces are over- 
whelmed and the system behaves as if there were no symme- 
tries. Hence, the non-chaotic QKT reverts to non-exponential 
fidelity decay. RMT predictions are also shown (dashed lines). 
All plots are averaged over 100 random initial states. 

always exponential. This tells us that the fidelity decay 
due to a quantum error is exponential independent of the 
chaoticity of the system even when applied after every 
map iteration. Here there is no concern that the system 
has random eigenvectors in the basis of the perturbation 
since the perturbation is diagonal (and in this basis the 
regular QKT eigenvectors are not random). Initial angu- 
lar momentum coherent states for the QKT and random 
states for the quantum Harper's map give similar results. 

Exponential fidelity decay due to quantum perturba- 
tions is not universal. We demonstrate this via unitary 
matrices of the interpolating ensembles [2^, which are 
intermediate between random, 6 = 1, and diagonal with 
Poissonian distributed eigenangles 6 = 0. Matrix statis- 
tics for these ensembles transition smoothly between the 
two limits [2^ . There is no known way of efficiently 
implementing such matrices on a quantum computer. For 
our purposes matrices with 6 < 1 provide models of non- 
chaotic operators that have no classical analog. Any per- 
turbation of these operators must be quantum as, a pri- 
ori, the operators have no classical parameters. Fig. [21 
demonstrates that the 6 1 fidelity decay for the per- 
turbation of Eq. 21 e = .3, is not exponential. 

An explanation for the above discrepancies is the pres- 
ence of symmetries. Non-chaotic systems containing 
symmetries, and thus invariant subspaces, exhibit chaos- 
like exponential fidelity decay, while non-chaotic systems 
without symmetries exhibit faster-than-exponential de- 
cay. This result is surprising. A priori we would expect 
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FIG. 2: (Color online) Fidelity decay for matrices of the in- 
terpolating ensemble 5 = 1 (•), .98 (o), .94 (□), .9 (+), .8 
(o), .7 (x) for the collective qubit z rotation with A'' = 128 
and e = .3. The fidelity for 5 < .7 matrices has the same 
behavior as those with (5 = .7. For (5 7^ 1 the fidelity decay is 
faster than exponential even though the perturbation is not 
tied to a classical parameter (the matrices have no classical 
analogs). As 5 ^ 1 the matrices become random and the 
decay rate approaches the RMT rate (dashed line). All plots 
average over 10 operators with 100 random initial states per 
operator. The lower inset shows the local density of states 
for an average of 50 5 = 1, .9, .8, .7 operators, compared 
to the Lorentzian expected for chaotic systems (solid line), 
and a Gaussian (dashed line). The Gaussian like LDOS for 
non-chaotic systems shows that the perturbation is localized. 
The upper inset shows the fidelity decay for e = .3, A/^ = 256 
block diagonal operators in which the two A'^/2 = 128 blocks 
are difi^erent interpolating ensemble matrices of the same J. 
Changing bases and having the perturbation break the sym- 
metry of the operator, brings the fidelity decay close to the 
RMT prediction (dashed line). As 5 is decreased, the fidelity 
decay becomes slightly slower than exponential. 

systems with less symmetries to be more chaos-like. Be- 
low we provide further numerical evidence demonstrating 
this behavior: one QKT symmetry sector and interpolat- 
ing ensemble matrices with added symmetries. 

Following Peres (J we block diagonalize an even J 
QKT and keep one block as the system. Using the per- 
turbation of Eq. ^ the fidelity decay for the block is 
faster-than-exponential for non-chaotic values of fc and 
exponential for chaotic values of fc, as shown in Fig. O 
This is in contrast to the exponential decay seen for all fc 
values when using the complete QKT. Thus, by remov- 
ing the symmetry, the system deviates from the expected 
RMT behavior. In addition, for chaotic values of fc, the 
exponential decay rate is not the one given by Eq. 12 
again in contrast to the full QKT where the decay rate 
was exact. 

We expect the fidelity decay of systems with sym- 
metries and extremely strong perturbations, such that 
the symmetries are overwhelmed, to revert to non- 
exponential. This occurs in the full QKT with e = .6, 
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FIG. 3: (Color online) Fidelity decay for the QKT subspace 
even with respect to 180° rotations about y and odd with re- 
spect to 180° rotations about x, J = 1024, A^ = 512, fc = 1 
(o), 3 (□), 5 (o), and 7 (A) for the collective bit z rotation 
with e = .2. Though the perturbation is not attached to a 
classical parameter the fidelity decay for non-chaotic QKTs 
are non-exponential. The chaotic QKT exponential decay 
does not follow the RMT rate, F = .45 (dashed line, actual 
rate is F ~ .405). The inset shows the fidelity decay for the 
same parameters with the perturbation transformed into a 
random basis. The decay is exponential for all values of fe at 
the RMT rate. 

Fig. n] inset. 

A symmetry can be added to interpolating ensemble 
matrices by using two matrices of equal b and Hilbert 
dimension N jl as diagonal blocks oi a. N x N operator. 
The perturbation of Eq. ^does not cause mixing between 
the blocks so we transform into a new basis. The first 
transformation we apply takes the perturbation to a col- 
lective qubit X rotation, ax replaces CTz in Eq. 0] With 
this the fidelity decay is exponential at the RMT rate 
even for low d blocks. However, this exponential decay is 
not due to the broken symmetry. Rather, the matrices of 
the interpolating ensembles, though not random with re- 
spect to the ffz collective qubit perturbation, are random 
in the eigenbasis of the collective qubit perturbation. 
As explained in 6] (and shown in Fig. O for the QKT 
without symmetry), when the system is random in the 
perturbation basis, the fidelity decay is exponential. 

A transformation which does not induce randomness 
in a block diagonal operator with interpolating ensemble 
blocks is a modified version of the transformation matrix 
to block diagonafize the even J QKT 0, [sOl ■ With this 
transformation the fidelity decay is exponential for S ^ 1 
blocks and becomes slightly slower than exponential for 
lower values of 5, Fig. |21 Thus, by adding a symmetry to 
the system we have slowed the fidelity decay and made 
its behavior more chaos-like. 

The above simulations demonstrate that perturbations 
of non-chaotic systems, even ones not attached to a clas- 
sical system parameter, tend to be local. However, when 
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the perturbation breaks a symmetry there are long range 
effects similar to those of chaotic systems. This process 
is different than that described in Ref. where the ex- 
ponential decay is linked to the RMT statistics of the 
system eigenvectors in the perturbation basis. 

Coherent far-field errors, such as Eq. 0] are the sub- 
ject of many theoretical and experimental quantum error- 
correction codes and encodings [27I Isif . Our results sug- 
gest that this error may cause only exponential fidelity 
decay even when a quantum computer is simulating non- 
chaotic evolution. This provides a novel way of protecting 
against these noise operators, add a symmetry to the sys- 
tem. Such symmetries already exist in encoded qubits, 
where specific states of a multi-qubit system are used as 
a logical qubit. For example, logical qubits in quantum 
dots that are encoded such that the exchange coupling is 
universal jsS, are symmetric with respect to the Sz 
angular momentum operator. In general, by enforcing all 
reasonable decoherence mechanisms to break a symme- 
try, decoherence may cause an exponential, rather than 
Gaussian decay. 

Exponential decay occurs in many physical systems not 
generally regarded as chaotic, including relaxation phe- 
nomena and Fermi Golden Rule calculations. We sug- 
gest, without proof, that symmetry may add insight into 
the preponderance of exponential decay laws in nature. 
Lack of symmetry may explain why some systems, such 
as magnetic relaxation in single molecular magnets |34|. 
decay non-exponentially. 

In conclusion, we have studied the effect of symme- 
tries on fidelity decay behavior. When the perturbation 
is not tied to a classical parameter of the system, as would 
likely arise in quantum computers, the presence or lack 
of symmetries strongly affects the fidelity decay behavior. 
Surprisingly, the presence of symmetries forces the sys- 
tem into a more chaos- like behavior, exponential decay, 
while lack of symmetries causes deviations from RMT 
predictions and a faster-than-exponcntial decay. Build- 
ing symmetries into quantum computers, as done when 
encoded qubits, can cause decoherent processes to affect 
the system with the less-damaging exponential decay. 
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